Comment on "Torsion Cosmology and the Accelerating Universe" 
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Cosmological solutions for homogeneous isotropic models in the framework of the Poincare gauge 
theory of gravity based on gravitational Lagrangian adopted in the paper by Kun-Feng Shie, James 
M. Nester and Hwei-Jang Yo (Phys. Rev. D 78, 023522 (2008)) are discussed. Cosmological 
solutions for accelerating Universe obtained in referred paper are compared with corresponding 
solutions of standard ACDM-model and also with cosmological solution obtained by authors of this 
Comment. 

PACS numbers: 04.50.+h; 98.80.Cq; 95.36.+X 



The referred paper [H is devoted to investigation of 
isotropic cosmology in the frame of the Poincare gauge 
theory of gravity (PGTG) based on gravitational La- 
grangian, which contains besides scalar curvature also 
the term quadratic in the scalar curvature and three 
terms quadratic in the torsion tensor with some coeffi- 
cients. The PGTG with such gravitational Lagrangian 
was adopted in previous works of authors as describing 
dynamical torsion field. According to [l| corresponding 
theory can give the solution of the problem of the accel- 
eration of cosmological expansion at present epoch (dark 
energy problem of general relativity theory (GR)). Note 
that this problem was investigated in Q in the frame 
of PGTG by using general expression for gravitational 
Lagrangian including both a scalar curvature and var- 
ious terms quadratic in the curvature and torsion ten- 
sors with indefinite coefficients, where regular inflation- 
ary Big Bang scenario with accelerating stage of cosmo- 
logical expansion at asymptotics was proposed. Because 
the behaviour of cosmological solutions for accelerating 
Universe obtained in [l| and are essentially different, 
we will analyze below in what relation to cosmological 
solution obtained in 0] is the cosmological solution of 
With this purpose we will consider the applica- 
tion of gravitational equations of PGTG for homogeneous 
isotropic models (HIM) deduced in in the case of grav- 
itational Lagrangian used in [l|. 

The following expression for gravitational Lagrangian 
of PGTG was used in (by applying notations and def- 
initions of 0) [ll|: 

^ foF + F"^^" (/i FaPfj,^ + /2 Faf_il3u + fs Ff^yafi) 

+F^'' (/4 F^. + h F,^) + /e F2 
+5"^'' (ai S^^, + 02 S,^^) + aa 5%c<5/3'^^.(l) 

In general case any HIM in PGTG is characterized by 
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three functions of time: the scale factor of Robertson- 
Walker metrics R and two torsion functions Si and 5*2. 
Then not vanishing components of the curvature tensor 
can be expressed by the following 4 functions: 
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A3 = 2 (if -2^1) 52, 
A4 = S2 + HS2, 



(2) 



where H = R/R is the Hubble parameter and a dot de- 
notes the differentiation with respect to time. The total 
system of gravitational equations of PGTG for HIM filled 
with spinless matter with energy density p and pressure 
p takes the following form 0: 
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where 



Ai + 2H {Ai - A3) + ASiA. 



Si = 0, 



4/52A2 



(5) 



-2(ji52Ai-(/o-6)52 = 0, (6) 
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a = 2ai + 02 + 3a3, 6 = 02 — ai, 

/ = A + y + /3 + /4 + ./.5 + 3/6 , 
9l = /2-2/3 + /4 + /5 + 6/6, 92 = 2/1-/2. 
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The system of gravitational equations (3)- (6) allows to 
obtain equations for the torsion functions and S2, and 
also the generalization of Friedmann cosmological equa- 
tions for HIM in the frame of PGTG in the case of differ- 
ent particular gravitational Lagrangians. By analyzing 
the system of equations (3)-(6) we will use the Bianchi 
identities in the Riemann-Cartan continuum, which are 
reduced in the case of HIM to two following relations [§| : 



A2 + 2H{A2- Ai)+ASiAi+2S2Ai^0, (7) [| 



reduced to three equations (3)- (5) including two indefi- 
nite parameters: a and /. Note that these equations take 
place always independently on restrictions for indefinite 
parameters of gravitational Lagrangian (1), if 5*2 = 0. 
At first time these gravitational equations for HIM were 
deduced and investigated in The analysis of gravi- 
tational equations (3)-(5) by using the Bianchi identity 
(7) (with S2 = 0) gives the following expressions for the 
curvature functions ^2, Ai and the torsion function Si 



+ 2i/ (A3 - A4) + 451A4 - 252^1 = 0. (8) -^=^(-^0 + ^" 



In the case of gravitational Lagrangian used in [l[ the 
following restrictions on indefinite parameters in equa- 
tions (3)-(6) arc valid: (72 = and 2f ~ qi. As result by 
using (8) we obtain from (6): 



[2fF + 3(/o - b)] S2 = 0, 



(9) 



where F = 6{Ai +A2) is the scalar curvature. According 
to (9) there are two possibilities corresponding to two 
different cosmological solutions: 



82^0 and F 
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In the first case, when conditions (10) are valid, the 
system of equations (3)- (5) by using (7) and the defini- 
tions (2) for the curvature functions A2 and Ai leads to 
the following solution: S*! = and 
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The expression (12) for the torsion function 5*2 and cos- 
mological equations (13)-(14) are identical to that at 
asymptotics obtained in These equations allow to 
explain observable acceleration of cosmological expan- 
sion by certain restrictions on indefinite parameters 
(see also Q); however, because the equations (13)-(14) 
have the form of Friedmann cosmological equations with 
cosmological constant, the problem of cosmological sin- 
gularity remains in such theory. 

The second possibility, when conditions (11) are valid, 
is considered in [l|. In this case gravitational equation 
(6) vanishes, and the system of gravitational equations is 



Si 



= -\ (/o + \ 
Ai = \F-A2 



(/o + la^O). 



(15) 

Unlike GR, where gravitational equations for HIM are 
identical to Friedmann cosmological equations, discussed 
gravitational equations of PGTG for HIM include besides 
the scale factor R also the torsion function Si. How- 
ever, gravitational equations allow to deduce cosmologi- 
cal equations for the scale factor R without torsion func- 
tion generalizing Friedmann cosmological equations. We 
obtain these equations by substituting the solution (15) 
into the definitions (2) of the curvature functions A2 and 
Ai. The explicit form of the first cosmological equation 
is the following: 
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(16) 



where the scalar curvature F is defined by (15). We 
do not write the second cosmological equation following 
from (16) by using the conservation law, which in the 
case of spinless matter minimally coupled with gravita- 
tion takes the same form as in GR: 



p + 3H{p + p) = 0. 



(17) 



Cosmological solutions for metric quantities can be found 
by solving the system of equations (16)-(17) by given 
equation of state for gravitating matter and initial con- 
ditions for (i?, H, H, p). At first note that if the pa- 
rameter a is not equal to zero (in [l[ a = Aa[ > [11]), 
cosmological equations contain higher derivatives of the 
scale factor i?; in particular, the cosmological equation 
(16) contains the third derivative of R (the second cos- 
mological equation is differential equation of the fourth 
order with respect to R) p^ . By given initial condi- 
tions we obtain two different solutions corresponding to 
two different initial values of H (or 6*1). It is because 
the equation (16) contains the squared derivative H and 
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hence by given initial values of {R, H, H, p) there is 
two different values of H . Cosmological equations con- 
tain the parameter \-, — —-, — ~ a = > with 

{fo + ia) J la 

inverse dimension of energy density, and they transform 
into Friedmann cosmological equations of GR if / 
(a — > 0) The solutions asymptotics of cosmological 
equations, where energy density is small ap <C 1, is close 
to the Friedmannien asymptotics . The behaviour of cos- 
mological solutions at present epoch depends essentially 
on the value of parameter a. If the value of cor- 
responds to the scale of extremely high energy densities 
and hence at present epoch we have ap -C 1, the be- 
haviour of cosmological solutions is quasi-Fricdmannicn 
and such solutions can not describe cosmological acceler- 
ation. If the value of a^^ is comparable with the average 
energy density at present epoch [ap 1), the solutions 
behaviour differs from that of GR and depends essentially 
on values of indefinite parameters (/, a) and using initial 
conditions. With the purpose to compare cosmological 
solutions for flat model (fc = 0) with that of standard 
ACDM-model of GR, we will use the initial values for 
{H, H, p) obtained in the frame of GR and correspond- 
ing to observational data. Note that as result by given 
parameters / and a the initial values for F and Si are de- 
termined and they can not be introduced independently 
from / and a (compare with [l|). 

Now let us to demonstrate these statements by nu- 
merical solution of Eqs. (16)-(17) in the case of flat 
model filled with the dust matter (p = 0), for which 
pR^ = const. With this purpose dimcnsionlcss units 
for variables and parameters will be introduced by the 

following way: t ^ t = iy^, H ^ H = H,J^, 
H H' = ^ p = -P-n a = ^ and 

dt ^ Per ' JO 

f ^ f ~ f-Tfi, where the value of pcr = 6/o-ffn cor- 
responds to average energy density in the Universe at 
present epoch in the frame of GR (the index " 0" at physi- 
cal quantities denotes its values at present epoch) . In the 
frame of standard ACDAf-model one supposes that the 
total value of energy density includes the contributions 
of three components: baryonic matter, dark matter and 
dark energy. By taking into account that one uses at 
present epoch for baryonic and dark matter the equation 
of state of dust (j>bm = Pdm = 0), and for dark energy 
Pde = — Pde, we write the Friedmann cosmological equa- 
tions of GR in dimensionless form for considering case: 

= PBM + PDM + PDE, 

i^' + i?' = -^(/5bm + Pdm-2pde), (18) 

In accordance with Eqs. (18) and observational data we 
have the following initial values for physical parameters: 
Hq = 1, pbmo + PUMO = 0.3, ii'a = -0.45. If one sup- 
poses that there is in the Universe only baryonic and dark 
matter (p = p-QM + Pdm), then we obtain the following 
initial condition for energy density po = 0.3. 

By using obtained initial conditions numerical solu- 
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FIG. 1: Solution for the Hubble parameter by initial condi- 
tions: po = 0.3, Ho = 1, Ho = -0.45, Hq = -4.60 and 2 = 2, 
/ = 0.35. 
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FIG. 2: Solution for acceleration w by initial conditions: pb = 
0.3, Ho = l, H'o = -0.45, H'o = -4.60 and a = 2, / = 0.35. 



tions for dimensionless Hubble parameter and accelera- 
tion w ^ H' + are presented in Fig. 1-2 for the 
following choose of parameters: / = 0.35 and a = 2. 

Similar solutions were used in [l|. Note that the solu- 
tion presented in Fig. 1-2 corresponds to the following 
initial value of H" (or ^1): i/^' = -4.60 (5*10 = 0.185). 
By using the second initial value of i^Q = —32.5 (S'lo = 
0.815) compatible with used initial values for i7, H' and 
p, we obtain not physical solution with non- Friedman- 
nien asymptotics - with rapidly growing values of energy 
density and the Hubble parameter in the future. Al- 
though by virtue of used initial conditions the solutions 
similar to that presented in Fig. 1-2 have at present 
epoch the same observational characteristics as solutions 
of the standard ACDM-model, however, the behaviour 
of cosmological solutions in the past as well as in the fu- 
ture differs essentially from that of GR. Unlike KCDM- 
model, where the transition from decelerating to acceler- 
ating expansion (w = 0) takes place at certain moment 
in the past, in the case of discussed solutions the accel- 
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FIG. 3: Solution for the Hubble parameter by initial condi- 
tions: po = 0.05, Ho = 1, H'o = -0.45, H'^ = -7.97 and 
a = 2, /= 0.35. 



w 




t 



FIG. 4: Solution for acceleration lu by initial conditions: po — 
0.05, Ho = 1, -ffo = -0.45, ifo = -7.97 and a = 2, / = 0.35. 

eration function w oscillates near w = 0. The character 
of solutions of eq. (16) with quasi- Friedmannien asymp- 
totics in the future depends essentially on initial values 
for H, H' and p. As illustration of this statement the so- 
lution of (16)-(17) corresponding to the following initial 
conditions: po = 0.05, Hq = 1, = -0.45, i?^' = -7.97 
(S'lo = 0.261) is presented in Fig. 3-4. Such initial condi- 
tions correspond to supposition that at present all matter 
in the Universe is reduced to the baryonic matter. We 
see that according to Fig. 3 the evolution in this case 
has the oscillating character. The second solution cor- 
responding to used initial values of p, H , H' is also not 
physical solution with non- Friedmannien asymptotics. 

We sec that in the case a ^ the character of cos- 
mological evolution differs essentially from that of GR. 
This is connected with the presence in cosmological equa- 
tions of higher derivatives of R. In order to exclude 
higher derivatives from cosmological equations the condi- 
tion a = was used in Refs. [1] and Then according 
to (15) the torsion 5*1 is certain function of energy density 
and pressure, and all cosmological solutions of general- 



ized cosmological Friedmann equations obtained in are 
physical solutions with usual Friedmannien asymptotics. 
These equations allow to build regular isotropic cosmol- 
ogy including inflationary cosmology 0, S 13 • Regular 
character of all cosmological solutions with respect to 
metrics, Hubble parameter, its time derivative and en- 
ergy density is connected with gravitational repulsion 
effect provoked by the torsion function Si at extreme 
conditions in the beginning of cosmological expansion. 
However, the torsion function Si tends to zero at asymp- 
totics, and such theory can not describe the accelerat- 
ing Universe. As was shown in regular HIM with 
quasi-Friedmannien accelerating cosmological expansion 
at present epoch contain two torsion functions 5*1 and 
S2, and the effect of gravitational repulsion at present 
epoch is connected with effective cosmological constant 
induced by the torsion function S2 ■ As it follows from our 
consideration given above, the gravitational Lagrangian 
adopted in Q allows to build HIM only with one torsion 
function. In the case of HIM with 5*2 = the scalar cur- 
vature F is not constant for physical solutions (see (15)). 
However, if the dominator in (15) is equal to zero (that 
corresponds to constant scalar curvature), gravitational 
equations do not allow to determine the torsion function 
Si and the curvature functions (compare with [l[ ) . Grav- 
itational equations have in this case the following specific 
solution: cosmological equations take the form of Fried- 
mann cosmological equations with cosmological constant, 
in which instead of Newton's gravitational constant fig- 
ures the parameter inversely proportional to (—a). If 
a > 0, this specific solution is not physical. In connec- 
tion with this the sign of a for this case was changed 
that leads to "degenerate" case discussed in [l[. Note 
that specific not physical solution with constant scalar 
curvature does not appear if a = 0. 

The investigation of isotropic cosmology in the frame 
of PGTG based on gravitational Lagrangian (1) shows 
that the PGTG can have principal meaning for theory 
of gravitational interaction Physical consequences 

depend essentially on restrictions on indefinite parame- 
ters of Cg. Because at present the quadratic part of Cg 
is unknown, wc have to investigate the PGTG by using 
general expression (1) of gravitational Lagrangian and to 
obtain restrictions on indefinite parameters of Cg in order 
that physical consequences and mathematical structure 
of gravitational equations should be satisfactory. 

As it is follows from our consideration given above the 
PGTG based on gravitational Lagrangian used in Q (and 
also in [T^]) lead to isotropic cosmology, in the frame 
of which the evolution of the Universe differs essentially 
from that of standard ACiI'M-model. Unlike [l[ cosmo- 
logical equations deduced in take at asymptotics the 
quasi-Fricdmannien form, and terms related to dark en- 
ergy (and possibly to dark matter Q) in cosmological 
equations of general relativity theory for ACDAf -model 
are connected in considered theory with the change of 
gravitational interaction provoked by spacetime torsion. 

One of the authors (A.V.M.) is very grateful to Profes- 
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